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Abstract. We present a systematic analysis of the description of odd nuclei by the Skyrme-Hartree-Fock 
approach augmented with pairing in BCS approximation and blocking of the odd nucleon. Current and spin 
densities in the Skyrme functional produce time-odd mean fields (TOMF) for odd nuclei. Their effect on 
basic properties (binding energies, odd-even staggering, separation energies and spectra) is investigated for 
the three Skyrme parameterizations SkI3, SLy6, and SV-bas. About 1300 spherical and axially-deformed 
odd nuclei with 16 < Z < 92 are considered. The calculations demonstrate that the TOMF effect is gener- 
ally small, although not fully negligible. The influence of the Skyrme parameterization and the consistency 
of the calculations are much more important. With a proper choice of the parameterization, a good de- 
scription of binding energies and their differences is obtained, comparable to that for even nuclei. The 
description of low-energy excitation spectra of odd nuclei is of varying quality depending on the nucleus. 



PACS. 21.60.Jz Nuclear density functional theory 
21.10.Pc Single-particle energies 



21.10.Dr Binding energies and nuclear masses 
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1 Introduction 

Nuclear density-functional-theory (DFT) in its most promi- 
nent versions, Skyrme-Hartree-Fock (SHF), Gogny, and 
relativistic mean- field (RMF), are widely used theoretical 
tools for the description of ground states and dynamics of 
atomic nuclei, see the reviews [TJOE]. They have reached 
a high level of quality in description of nuclear properties 
throughout the nuclear chart and, what is important, do 
this in a self-consistent manner. The overwhelming major- 
ity of applications deal with even-even nuclei which are far 
simpler to handle than odd ones. Odd nuclei add crucial 
pieces of new information, e.g. single-particle spectra, and 
thus are increasingly studied in recent years. A fully self- 
consistent calculation for odd nuclei has to account for 
all the terms following from the initial functional. This 
means, in particular, that contributions of both time-even 
and time-odd densities from the initial functional are to 
be used [1,2,3,4,5 . The time-odd densities (current j, 
spin s, and vector kinetic-energy T) and related time-odd 
mean fields (TOMF) do not contribute to the ground state 
of non-rotating even-even nuclei but can be essential for 
the ground states of odd and odd-odd nuclei, as well as 
in nuclear dynamics including nuclear rotation 01(5,7,8 
and electric [ MTOlfTTl[l"2"] . and magnetic [TBlfTl] giant res- 
onances. Moreover, TOMF embrace the crucial terms to 
restore the local gauge invariance (Galilean invariance in 
non-relativistic Skyrme and Gogny models) violated by 
velocity-dependent time-even densities |H[S]. 



In odd nuclei, TOMF contribute to the ground state 
mean fields and, by breaking the time-reversal symme- 
try in the intrinsic frame, destroy the Kramer's degener- 
acy of the single-particle states. TOMF might also influ- 
ence the single-particle spectra and such characteristics 
as binding energies, odd-even staggering, and separation 
energies. Note that binding energies and separation en- 
ergies are important for astrophysical applications [3lll5|. 
the odd-even staggering determines the quality of the pair- 
ing description [16U17IITgl lT! P P T P 2 ] . the separation en- 
ergies are crucial for estimation of nuclear stability in drip 
line [55] and super-heavy regions [23] , and the low-energy 
single-particle spectra are now in focus of exploration of 
exotic nuclei [2"3"I|2"5] . 

The role of TOMF in odd nuclei has been already 
thoroughly inspected for RMF models, see e.g. the re- 
cent study [35] and references therein. Most applications 
of SHF models neglect TOMF. There are a few early pub- 
lications 18,27 taking care of TOMF, however, involving 
only light nuclei and a limited number of Skyrme forces 
(SIII [28 and SLy4 [2H])- The large interest on a proper de- 
scription of odd nuclei motivates a fresh look at the case. 
Thus very recently, a new SHF study of rare-earth odd 
nuclei has been published [3D] , where the TOMF effect on 
excitation spectra of odd nuclei was discussed in detail. 
Altogether the previous studies have found rather weak 
influence of TOMF on properties of odd nuclei. The max- 
imal effect was obtained in light nuclei 1181126] . For bind- 
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ing energy, the TOMF were shown to give more binding 
within RMF [26] and effects of both sign within SHF [18]. 
The TOMF-induced changes in the single-particle spectra 
were found very small |26ll3"U] . It was claimed that RMF 
results much less depend on the force parameterization 
than SHF ones [2"6] , 

The previous studies were limited by particular mass 
regions or specific features of odd nuclei. At the same, it 
would be very instructive to have a general view of the 
problem, covering most of the nuclear chart and involving 
all the basic characteristics which might be affected by 
TOMF. Just this aim is pursued in the present paper de- 
voted to a systematic exploration of TOMF effects in SHF 
calculations for proton and neutron odd nuclei. The study 
covers about 1282 odd nuclei with 16 < Z < 92 (from 
sulfur to uranium) where experimental energy data are 
available. Both spherical and axially deform nuclei are in- 
volved. Three Skyrme parameterizations SkI3 [31] . SLy6 
[25] , and SV-bas [35] with essentially different effective 
masses (m*/m=0.58, 0.69, and 0.9, respectively) are ap- 
plied. A variety of observables (binding energies, odd-even 
staggering of energies, separation energies, and low lying 
single-particle spectra) are analyzed. To separate the ef- 
fects from time-odd spatial current j and spin density s, 
the time-odd terms in the Skyrme functional are sorted 
into the minimal set for restoring Galilean invariance, and 
the sets with additional spin and spin-gradient couplings. 
The results are compared with previous studies in SHF 
QI1[3Q] and RMF [26]. 

The paper is outlined as follows. In Sec. [21 the the- 
oretical framework and calculation scheme are outlined. 
In particular, the Skyrme functional with the time-odd 
terms and the pairing functional are presented. Results 
of the calculations are discussed in Sec. [3J The summary 
is given in Sec. 3] Some details concerning the functional 
parameters and single-particle wave functions are done in 
Appendices A and B. 



2 Model and details of calculations 



The starting point is the total energy E of system [T][3] 
E = f dr{Uki n (T q ) +H S k(pq,T q ,s q ,j g ,J q ,T q ) (1) 

+ 'Hpair(Xg) + Uc{Pp)} ~ E cm , 

involving local time-even (nucleon p q , kinetic-energy r q , 
spin-orbit J q ) and time-odd (current j q , spin s q , spin kinetic- 
energy T q ) densities, as well as the pairing density \q- The 
label q denotes protons and neutrons. 



Total densities, like p — p p + p n , have no the index. 
Further 

^ = £^(,), (2) 

9 H 
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are kinetic-energy, pairing, Coulomb, and center-of-mass 
terms. Note that the center-of-mass energy E cm is not 
included in the mean-field equations but evaluated a pos- 
teriori (32) . The pairing interaction may be density depen- 
dent (surface pairing with F = 1 — p/po and po = 0.20113 
fm -3 ) or not (volume pairing with F = 1) [T9"] . 

The key part of ([T]) is the Skyrme functional Hsk com- 
posed from the terms 
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where bi, b[, bi, b[ are the force parameters. Their con- 
nection with the standard Skyrme parameters ti and Xi is 
specified in the Appendix |XJ The dominant Skyrme con- 
tribution "H^ en ^ involves only time-even densities p q , r q , 
and J q . It embraces all standard terms relevant for the 
ground state properties and electric excitations of even- 
even nuclei (for electric modes the current can be also 
important 9. 10HT2"]). The isovector spin-orbit interaction 
is usually linked to the isoscalar one by b' 4 = 64 and the 
tensor spin-orbit terms oc b\, b[ are often omitted [1,3 . For 
the sake of simplicity we also omit them and the present 
study and so use b\ — and b[ — throughout. 



The next Skyrme contributions, W 



(Gal) ^(s 2 ) ^(sAs) 



Sk 



« sk , ft 



Sk 



add gradually terms with the time-odd densities j ? , s q , 
and T q and create corresponding TOMF in the mean- 
field equations. These are the terms of our particular in- 
terest. They do not contribute to ground state properties 
of even-even nuclei but become relevant in odd nuclei and 
nuclear dynamics (rotation, electric and magnetic giant 
resonances). In the present study we consider three op- 
tions for the TOMF impact: 
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The first option deals only with Hg^^ > i- e - involves only 
time-even densities. This violates Galilean invariance in 
nuclear dynamics [3] . Adding the term restores this 

invariance. This yields "Hg™ m \ the minimal set to achieve a 
consistent dynamical model. Derivation of the functional 
from a Skyrme force [531154] also yields the spin terms 

and Hg^ As ' ) with the parameters bi and b[ uniquely 
related to the bi, b[, see Appendix IA1 Following our expe- 
rience, the term leads to unstable ground states for 
many medium and heavy nuclei. This is a common prob- 
lem in most Skyrme parameterizations [551136] . Thus we 

include to the option % s ™ m+s ' only the simple spin term 



Hg k . The comparison of results for Hg^ LL) with those 

from Hg™ m ^ and % s ™ m+s ' should provide information on 
the influence of TOMF. The corresponding total energies 
are denoted as E even , E min , and E min+S 2. 

The calculations are performed in cylindrical coordi- 
nate space grid with the mesh size 1 fm. The axial equi- 
librium deformations are determined by minimization of 
the total energy [1]. The calculations are restarted from 
different initial deformations to separate isomeric and ab- 
solute minima of the energy. 

To explore the possible variance of the results, the cal- 
culations are done for three different Skyrme parameteri- 
zations, SkI3 [31], SLy6 [55], and SV-bas [3J with the ef- 
fective masses 77i*/m=0.58, 0.69, and 0.9, respectively. All 
the parameterizations provide a reasonable description of 
stable even-even nuclei and giant resonances, though per- 
form differently with respect to more detailed observables. 
The force SkI3 is known to reproduce peculiarities of Pb 



(even) 



isotopic chain. SLy6 was developed with special empha- 
sis on neutron rich nuclei and neutron matter. The recent 
parameterization SV-bas [53] is fitted for a large data set 
covering long isotopic and isotonic chains and so may be 
most suited for systematic explorations. Note that SkI3 
and SV-bas allow an isovector spin-orbit force (decoupled 
parameter 6 4 ) while SLy6 fixes that to the isoscalar term 
by setting b' 4 = 64. Both these differences can affect the 
single-particle shell structure and so the properties of odd 
nuclei. 

The pairing functional ([3]) uses the volume pairing for 
SkI3 and SLy6 and surface pairing for SV-bas. The pairing 
strengths \-p a ir,g are adjusted following the recipe [32]. The 
pairing is augmented by the stabilization procedure from 
37J. This prevents the breakdown of the pairing near the 
closed shells which to some extent simulates the similarly 
smoothing in particle-number projection and so avoids un- 
physical kinks in properties along collective deformation 
paths. 

The pairing mean-field equations are derived from the 
given pairing functional and solved in the BCS approxi- 
mation. The pairing density in ([3]) is 



V/( r ) = ~/2 u kVk 



keq 



er=± 



(13) 



where the sum runs over the single-particle states k, the 
state k is time-conjugate to k, Vk and Uk are Bogoliubov 
amplitudes, <pf^{r) is the spin cr-component of the single- 
particle wave function in cylindrical coordinates (see ap- 
pendix |B]). 

The weight in flT5J 



f'-i 



1 + exp {[0.5(e fc + e s ) - A, - AE q ]/n q } 



(14) 



is an energy- dependent cut-off factor limiting the pair- 
ing space to the vicinity of the Fermi energy. The cut-off 
is necessary to suppress unphysically large contributions 
from high-lying states, caused by the zero-range form of 
the pairing interaction 1 . We use here the recipe [T§] with 
the cut-off parameters amounting in average to the typical 
values AE q =5 MeV and /i 9 =0.5 MeV [35]. The value e fc is 
the single-particle energy of state k and X q is the chemical 
potential. 

In odd nuclei, the odd nucleon misses a partner to join 
the pairing scheme and thus blocks the occupied single- 
particle state. The blocking results in the one-quasiparticle 
state [3Jin] 



£|BCS) 



H (u k +v k a+a+)\0) , (15) 



fe>0,fc#fci 



where a, 



and at are quasiparticle and single-particle cre- 



ation operators; |BCS) and |0) are BCS and HF vacuum 
states. The blocking effect is known to be of a crucial im- 
portance for odd nuclei [39ll40ll4T] . 

We compute the blocked states for a variety of single- 
particle states k near the Fermi energy. The state with 
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lowest energy is treated as the ground state while oth- 
ers form the low-energy excitation spectrum of the odd 
nucleus. 

Finally note that time-conjugate states are easily ob- 
tained in even-even nuclei as <j>^ = Tipk by using the time- 
reversal operator T. However, this is not trivial for odd nu- 
clei where the time-reversal invariance is broken. In this 
case, <j>j: is found as the state having maximum overlap 

with Ttfik- 



3 Results and discussion 

In this section we present our results for the influence of 
TOMF on binding energies, odd-even staggering, separa- 
tion energies, and low-energy spectra of odd nuclei. Three 
options (fllj]) - (ll"2l with the corresponding total ground state 
energies, E cvcni E min , and -E m i n +s 2 j are compared. The 
same index convention is applied to other observables. 

In this study we neglect the coupling with vibrational 
states of the even-even core [4TJHT] and nuclear rotation 
[42,43]. The loss of the corresponding collective correla- 
tions can result in underestimation of the total binding 
energies by 1-2 MeV in soft or deformed nuclei |44ll4"5] . 
For the variables represented through the differences of 
the energies (pairing gaps and separation energies) the ef- 
fect should be much smaller. For single-particle energies 
the role of the vibrational and rotational coupling can be 
dramatic [12] • In general it is reduced to much more com- 
pressed low-energy spectrum of odd nuclei |40U41] . So we 
may expect from very beginning for more dilute SHF spec- 
tra as compared to the experimental data. 



3.1 Binding energies 

Results for the binding energies (BE) are given in Figs. 
Q3H] and Table 1. The experimental data are taken from 
[46] . To amplify the effects, the energy differences AE are 
considered: -E CVO n — E m i n to check the effect of restoration 
of Galilean invariance, E QVCn — E min+S 2 to demonstrate 
the total TOMF effect when spin terms s 2 are also added, 



and E, 



min+s^ 



E exp to estimate the deviation from the 



experimental data. 

Fig. [T] shows the BE differences for odd Sn isotopes. 
Following panel a), the restoration of Galilean invariance 
enhances the binding. The effect is small throughout and 
practically negligible for SV-bas. This is not surprising as 
SV-bas with its effective mass m* /m — 0.9 has smaller co- 
efficients &i and b[ than SLy6 and SkI3 and hence weaker 
current contribution ~ j 2 . The additional inclusion of spin 
terms (panel b) counterweights the extra binding and leads 
to its modest reduction by 0.1-0.3 MeV. Thus the effects 
of TOMF on the BE are generally small. 



Panel c) of Fig. 1 shows the difference E. 



min+s- 



E, 



cxp 



between theory and experiment in both even and odd Sn 
isotopes. This difference is much larger than the previ- 
ous ones and shows remarkable isotopic trends. The typ- 
ical patterns for SkI3 and SLy6 show a tendency to over- 
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Fig. 1. The energy differences E cvcn - E min (a), E cvcn - 
i? min+s 2 (b), and -E m i n+S 2 — E cxp (c), in odd Sn isotopes for 
the parameterizations SkI3, SLy6, and SV-bas. 



binding at the shell closures N=50 and 82 and under- 
binding at mid shell. The most recent parameterization 
SV-bas produces the smoothest trends and smallest devi- 
ations extending to the limits of known isotopes. This is 
due including the long isotopic and isotonic chains into the 
SV-bas fit |32) . Even here there remains a region of under- 
binding by ~ 2 MeV near N=55-65. But these nuclei are 
soft and, as was mentioned above, should additionally ac- 
quire just this amount of the collective correlation energy 
|32U44| . Anyway, we see large deviation from the experi- 
ment and strong dependence of the results on the Skyrme 
parameterization. Both these factors are an order of mag- 
nitude stronger than the TOMF effects shown in panels 
a) and b). 

Fig.[5]shows the BE differences for the chain of isotones 
with the neutron number N=82. Most of the patterns are 
much the same as for the isotopic chain in Fig. [T] and 
so the above remarks can be applied for Fig. [2] as well. 
The only exception seems to be a different behavior of 
E m i n + S 2 — E cxp for SkI3 and SLy6 in panels c) of Figs. 
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Fig. 2. The same as Fig. [T] but for N=82 isotones. 



1 and 2, where we see similar SkI3 and SLy6 trends for 
isotopes and opposite for isotones. Perhaps, this is related 
with different fit strategies of these two parameterizations. 

In the following we discuss nuclei all over the chart of 
isotopes, even-even nuclei and proton-odd or neutron-odd 
ones called henceforth "odd". Doubly odd nuclei are not 
considered. 

Fig. [3] shows the difference -E^vcn — £min+s 2 all over 
the nuclear chart with 16 < Z < 92. The effect of TOMF 
is weak everywhere, never exceeding 1 McV and usually 
remaining even much smaller. In general the TOMF lead 
to less binding. The effect is stronger in light nuclei and 
decreases with the nuclear size, which is in accordance 
with RMF findings [33] . There is some dependence on the 
Skyrme parameterization concerning the increase of the 
TOMF effect for small nuclei, which is more pronounced 
for SkI3 and less for SV-bas. Note that, in contrast to the 
previous Skyrme study [18] for light nuclei, neither the 
RMF [26] nor our calculations find an enhancement of the 
TOMF effect at N=Z. Perhaps in our case this is caused 
by omitting the term s • As which was taken into account 
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Fig. 3. The binding energy difference -E CVC n — -E-min+s 2 for the 
Skyrme parameterizations SkI3, SLy6 and SV-bas, drawn in 
the N-Z plane. All available neutron and proton odd nuclei 
with the charge 16 < Z < 92 are included. 



in [T8] and gave there a strong contribution opposite to 
one from s 2 -terms. 

Fig. 0] shows the systematics of the deviation of cal- 
culated BE from the experimental data |46] in even and 
odd nuclei altogether. As all three parameterization show 
qualitatively the same pattern, we show only the results 
from SV-bas. As might be expected, the best agreement is 
found along the isotopic and isotonic chains of semi-magic 
nuclei. These are the nuclei with least correlation effects 

The present calculations omit these effects and gen- 
erally yield under-binding for deformed nuclei. Subtract- 
ing proper correlation energies cures the problem for light 
and medium nuclei. But there remains a systematic trend 
to under-binding of deformed heavy and super-heavy ele- 
ments which is a well known problem of all Skyrme param- 
eterizations [32p47] . The same trends exist for even-even 
nuclei and for other Skyrme parameterizations. Actually 
they are features of the present SHF functionals as such. 
Most important for our purposes is that odd nuclei do not 
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Fig. 4. The difference between SHF and experimental values 
46 of the binding energy, -E mm+S 2 — -E exp , calculated for the 
force Skl3 and plotted in the N-Z plane. All even and odd 
nuclei with 16 < Z < 92 are included. 

Table 1. R.m.s. deviations (in MeV) between SHF and exper- 
imental data |46| for the binding energies in odd, even, and all 
(odd + even) nuclei from Figs. [3] and [4] 



Parameterization 


Option 




Nuclei 








even 


odd 


all 




E C vcn 


3.51 


3.57 


3.55 


SkI3 


Em in 


3.51 


3.50 


3.50 




E m in+s 2 


3.51 


3.73 


3.66 




E C vcii 


4.18 


4.01 


4.07 


SLy6 


Emin 


4.18 


3.97 


4.04 




E m in+s 2 


4.18 


4.14 


4.15 




E C vcn 


2.81 


2.72 


2.75 


SV-bas 


E m i n 


2.81 


2.71 


2.74 




E • , 2 
^mm + s z 


2.81 


2.84 


2.83 



change the picture at all. They would merge smoothly 
into the data from even nuclei without showing any spe- 
cial even-odd effect. Proper tuning of BE in even nuclei 
thus automatically provides good results for BE in odd 
nuclei as well. 

Additional information on the global quality of BE 
description is done in Table Q] in terms of r.m.s. devia- 
tion from experimental data |46j . The numbers have to 
be taken with care because they are computed from raw 
mean-field energies without rotational projection and vi- 
brational correlation energy. 

Table [1] shows that restoration of the Galilean invari- 
ance (step from E even to -Emin) improves the description 
only a bit. But adding the s 2 -terms, i.e. stepping to E min+S 2 , 
again enhances the deviations. Note that even and odd nu- 
clei are described about with a similar accuracy. SV-bas 
yields to a better description than SkI3 and SLy6 and, 
again, the difference between the Skyrme parameteriza- 
tions is much larger than the TOME effect. 

Finally it is worth to compare our results with the 
RMF predictions, in particular with the recent study of 
[2"B"] covering light nuclei with 18<Z<27as well as Ce 
and Cf isotopes. Like in our case, they predict a weakening 
of the time-odd effects with increasing nuclear size. At the 
same time, the RMF calculations in light nuclei demon- 
strate a profound TOMF effect near the proton drip-line 
and generally more binding induced by the odd-time fields. 
These findings are not confirmed by our calculations, in 



particular, the Skyrme calculation usually reduce binding 
when including spin terms. However, we should not overes- 
timate the significance of this discrepancy. In both RMF 
and Skyrme calculations the TOMF contribution to the 
binding energies consists of a terms of a different sign and 
actually the final result is determined by a subtle balance 
of these terms. Using another parameterization or adding 
still neglected contributions (like s • As in the SHF) may 
change the balance and so magnitude and sign of a week 
TOMF effect. 

In agreement to the previous SHF study [IB] , the present 
calculations demonstrate a strong dependence of the re- 
sults on the parameterization. This is to a large extent 
a consequence of using essentially different parameteriza- 
tions SkI3, SLy6, and SV-bas (which possess much dif- 
ferent effective masses and vary also in other features). 
The RMF results on the other hand, are claimed to 
be almost independent of the parameterization. But it is 
ought to mention that the RMF parameterizations in this 
study are rather similar (all from group A in the classifica- 
tion of [26]). More differences may appear when utilizing 
other RMF parameterizations which have higher effective 
masses. 

3.2 Odd-even staggering (OES) 

During last years the odd-even staggering (OES) was a 
subject of intense investigations within various models in- 
volving RMF [17l[26], SHF [TUHSUTUEniEIlES] , and Gogny 
[4"8Il4"9"] approaches. Being mainly attributed to the pair- 
ing, the OES is widely used as a measure of the pair- 
ing gap. However, the OES can be also affected by other 
physical mechanisms (e.g. nuclear deformation or TOMF) 
which significantly complicates the problem and calls for 
discrimination of the various influences. For SHF, the TOMF 
effect on OES was so far explored either indirectly [2UJ or 
only for light nuclei [TH]. To the best of our knowledge, 
the present study is the first direct and systematic inves- 
tigation of the TOMF influence on OES. The experimental 
values of OES used in the present exploration are obtained 
by using evaluations [46] , 

There exist several options for characterizing the OES 
[39] . We adopt the three-point formulas 

Af{Z,N) = t^_ [E {Z,N-l)-2E{Z,N) 

+ E(Z,N+1)] , (16) 

Af (Z, N) = t^l [ E {Z — 1,N) — 2E(Z, N) 

+ E(Z + l,N)}. (17) 
It ought to be mentioned that background contributions 

(3) 

act differently for centering A z N at even or odd nuclei. 
When centering at odd-N or odd-Z nucleus, respectively. 
The three-point formulae have the advantage of canceling 
the smooth background of mean field contributions. In- 
stead, when being centered on even nuclei, (fT6]) - ([T7]) are 
likely to contain also significant shell and deformation ef- 
fects [HEPJ. 
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Fig. 5. Neutron gaps A$ in Sn isotopes, calculated for the 
parameterizations SkI3 (panel a), SLy6 (panel b), and SV-bas 
(panel c) for the options (|10[) - (112[) as indicated and compared 
with the experiment [46] . 



Fig. [5] shows the neutron gaps along the chain of Sn 
isotopes. All three parameterizations, SkI3, SLy6, and SV- 
bas, quite accurately describe the OES in sizes and trends. 
The best agreement is seen in the mid shell regions where 
pairing dominates in forming the gaps and fit of the pair- 
ing strength is accomplished [32]. Larger deviations and 
a strong dependence on the parameterization are found 
at and around the shell closures (N=50 and 82). That is 
a point where shell effects dominate and the shell struc- 
ture varies substantially amongst the parameterizations 
(e.g. due to the much different effective mass). At a more 
detailed level, one can see that, for SkI3 and SLy6, E m i n 
slightly underestimates Aft' while E min+S 2 gives the oppo- 
site effect. For SV-bas, both E min and E min+S 2 results are 
close to experiment. The description is not so perfect for 
the N=82 isotones shown in Fig. |5]where all three SHF pa- 
rameterizations show somewhat larger deviations from the 
experimental data, especially E min+S 2 for SkI3 and SLy6. 
But even these deviations are yet quite moderate and one 
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Fig. 6. The same as in Fig. [5]for the proton gaps A% (Z,N) 
in N=82 isotones. 



may speak on a generally good description of OES. The 

trends like a general decrease of A N with N and increase 

of A^p with Z are reproduced equally well as in the study 
|21j with the special isospin pairing interaction. 

Figs. [7] and [8] show the systematics of the deviations 
between calculated A$ (Fig. [7J or A ( f (Fig. [SJ and the 
experimental data. The majority of nuclei stays in the 
"green" window having perfect agreement with the data. 
We see strings of larger deviations along some semi-magic 
chains and regions of well deformed and transitional nu- 
clei. These nuclei are known to have large collective corre- 
lations [33] and are likely to need them for improvement 
of the gap description. The size of the deviations differs 
very much amongst the three parameterizations and SV- 
bas again gives the best results. This is most probably due 
to the large effective mass (m*/m — 0.9) which seems to 
produce a more correct level density. 

The r.m.s. deviations for OES are summarized in Ta- 
ble [2] Comparison with Table Q] shows that the average 
deviation for OES is much smaller than for the binding 
energies, which reflects the fact that OES are less sensi- 
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Fig. 7. The difference between the theoretical and experimen- Fig. 8. The same as in Fig. [7J for the proton pairing gaps 



I (3) 



tal values [46] of the neutron pairing gaps, A K N>min+s 2 -^jv^pi 
for all available nuclei in the window 16 < Z < 92 and for the 
Skyrme parameterizations SkI3, SLy6, and SV-bas. 



tive to nuclear bulk properties and related errors. This is 
one more example for the fact that energy differences are 
often more reliably described than a total energy as such 
[52] . Like for the binding energies in Table [II the TOMF 
effect on OES is smaller than variations between the dif- 
ferent parameterizations. 



3.3 Relation between OES and spectral gap 



The OES (fTl)|) - (fTT)) are closely related to the pairing gaps 
which appear in the two-quasiparticle excitation spectra 
of the BCS or HFB calculations [35]. Such spectral gap 
can be denned as the average [15] 



A q = 



Efee 9 \u k v k A k \ 



Efceg \ u k y k\ 



(18) 



We use for this purpose A k = J dr(p* k {r)A q 4> k {v) as the 
state-dependent gap in even and odd nuclei. The weight 



l(3), 



Z,N). 



Table 2. R.m.s. deviations (in MeV) between SHF and ex- 



perimental data [46] for A N z from Fig. [7] and A z from Fig. 

E 





SkI3 


SLy6 


SV-bas 


iJ cvcn 


0.98 


0.63 


0.36 


Emin 


0.97 


0.67 


0.37 


Emin+s 2 


1.01 


0.60 


0.36 




SkI3 


SLy6 


SV-bas 


Ecvcn 
Emin 
Emin+s 2 


0.79 
0.82 
0.82 


0.69 
0.73 
0.66 


0.39 
0.40 
0.35 



u a v a concentrates in the pairing-active region near the 
Fermi energy [15] . It is the spectral gap in even nuclei 
which is used to calibrate the pairing strength by compar- 
ison with the experimental A q 3 ^ [32]. A fully consistent 
computation of odd nuclei allows to check the resemblance 



of the spectral gap and A 



(3) 
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Fig. 9. Neutron gaps in Sn isotopes for three different Skyrme 
forces. Compared are the experimental gap from OES A 



the theoretical gap from OES A N min+s2 , and the spectral gap 
as defined in eq. ()18p . 



N \exp ' 



Fig. [9] shows the spectral neutron gap An together 
with the theoretical OES A^- min+s 2 an d experimental A^ 
in Sn isotopes. There is an acceptable agreement between 
A n and A N min+s 2 in the mid-shell region but a large de- 
viation in trend and magnitude near shell closures where 
pairing shrinks and shell effects dominate. The magnitude 
of deviation depends on the parameterization, which con- 
firms its shell origin. The deviations in the mid-shell re- 
gion stay in a range of about 10% and justify the resem- 

blance of An and A y N min+g2 as a first guess. The theo- 
retical results reproduce nicely the experimental data in 
the pairing-dominated mid-shell regions, which is not sur- 
prising since just these regions were used for the fit of the 
pairing strength. 

Fig. [TU] gives the similar comparison for the proton 
gaps in N=82 isotones. There is again a huge mismatch 



between Az and A 



(3) 



near the Z—50 shell closure. 



However the agreement in the mid shell region is not as 
good as for the neutron gaps and strongly depends on the 
Skyrme parameterization. It is thus likely that shell effects 




50 55 60 65 70 



Fig. 10. The same as figure [9] for proton gaps in the isotonic 
chain N=82. 



continue to contribute even far off the shell closure. Note 
that the deviation is smallest for SV-bas. 

It is interesting to note that one sees clearly that the 
spectral gap A was fitted to the experimental OES in even 
nuclei which shows still good enough agreement even in 
this case having generally large deviations. 

The sensitivity of OES to the pairing strength is tested 
in figure [TT] The result corroborates the finding of the pre- 
vious two figures. The pairing does not affect the regions 
near shell closures but dominates mid shell. The changes 
in the gaps are rather constant in mid shell and exceed 
those in the pairing strengths (30% versus 20%). A simi- 
lar amplification was previously found in |50] . 



3.4 Separation energies 

The neutron and proton separation energies (SE) 
S N (Z, N) = E(Z, N) - E{Z, N — 1) , 



Sz(Z, N) = E(Z, N) - E{Z - 1, N) 



(19) 
(20) 



represent an important nuclear characteristics. They are 
sensitive to all aspects of the nuclear mean field and pair- 
ing. Since SE deal with odd or/and odd-odd nuclei, they 
can be influenced by TOMF. 
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in the isotones with N=82. Lower: A 



in the Sn isotopes. 



(3) 



Figs. [T21 and IT51 show the systematics of the deviations 
of SHF neutron and proton SE from the experimental data 
|46j . There are pronounced errors along isotonic chains for 
Sn and isotopic chains for Sz . They sensitively depend on 
the parameterization and are minimal for SV-bas. Obvi- 
ously this is a shell effect. SE show a large jump at shell 
closures and size of the jump is closely related to the OES 
magnitude near the closure. The OES in turn, depends 
sensitively on the effective mass of the underlying model. 
The differences near the shell closures in Figs. [T^] and [T3] 
are thus related to the effective masses of SHF parameteri- 
zations as well. The good performance of SV-bas indicates 
once more that its effective mass of 0.9 is favorable. The 
regions of well deformed nuclei (remote from the magic Z 
and N lines) show generally smaller deviations for all pa- 
rameterizations. Instead the transitional regions (close to 
the magic Z and N lines) exhibit larger errors. The latter 
is probably caused by missing the correlations which are 
known to be most strong for transitional nuclei. 

Table [3] summarizes the r.m.s. deviations for SE. A 
large difference in performance for the three parameteri- 
zations is seen, where SV-bas shows the best results. The 
difference between various options of the energy functional 
is smaller. In particular, E m j n+S 2 gives a slight improve- 
ment of the overall quality. Figs. [T2TfT51 and Table [3] alto- 
gether indicate that the choice of the proper Skyrme pa- 
rameterization is crucial to achieve a good quality of the 
SE description in all regions of the nuclear chart. Transi- 
tional nuclei may require the correlation effects. 
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Fig. 12. The systematics of difference between SHF and ex- 
perimental 46 neutron separation energies, "S^min+s 2 *S*N,exp, 
drawn in the iV-Z-plane for the parameterizations SkI3, SLy6, 
and SV-bas. 

Table 3. R.m.s. deviations (in MeV) between SHF and ex- 
perimental data for separation energies Sn,z from Figs. 

mm 



Sn 


SkI3 


SLy6 


SV-bas 


Eeven 


1.15 


0.83 


0.50 


Fl ■ 

J-Tnin 


1.16 


0.85 


0.50 




0.80 


0.60 


0.50 


S z 


SkI3 


SLy6 


SV-bas 


J-'cvcn 


0.96 


0.83 


0.53 


F- ■ 


0.98 


0.87 


0.54 


E m in+ S 2 


0.97 


0.82 


0.50 



3.5 Excitation spectra of odd nuclei 

The excitation spectrum of odd nuclei is closely related to 
the single-particle spectrum of the even-even neighbor. In 
a simple picture one even may deduce the single-particle 
spectra of the even-even nucleus from the excitations of 
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the nearby odd systems. We will now check in more detail 
the effect of TOMF on the excitation energies of odd nuclei 
and the performance of the simple estimate in terms of the 
single-particle spectra of the even-even reference system. 
The comparison with the experimental data will be also 
done. Here, however, one has to keep in mind that the 
low-energy excitation spectra of odd nuclei may be heavily 
influenced by admixture of vibrational states of the even- 
even core, resulting in in a general essential compression 
of the low energy spectrum. [401I41] . 

Figure [14] shows the excitation spectra in 207 Pb. This 
nucleus has a neutron hole with respect to 208 Pb. The 
simple excitation model is then deduced from the neutron 
one-hole (lh) and one-particle-two-hole (lp2h) excitation 
spectra in 208 Pb, as indicated in the figure. In this model, 
the low-energy spectrum is dominated by lh excitations 
but at higher energies the lp2h excitations become ener- 
getically competitive. Such prescription is widely used in 
various approaches, see e.g. [39ll40ll4T| . and it obviously 
does not include the TOMF effects. 

These effects are illustrated in Fig. [TJ] for the options 
even, mean, and min + s 2 where the TOMF are added 



Fig. 14. Excitation neutron spectra in 207 Pb for the parame- 
terizations SkI3 (a), SLy6 (b), and SV-bas (c). In every panel, 
the experimental data [51] are compared with the SHF results 
for the options (|10p - (ll2|l . Also the simple one-hole (lh) and 
one-particle- two-hole (lp2h) estimates from the single-neutron 
spectrum of 208 Pb are exhibited. The dotted horizontal line 
marks energy of the lowest collective state 3~ in 208 Pb. For 
better view, results for one and the same level are connected 
by dash lines. See text for more details. 



step by step. Unlike the simple model discussed above, 
these results are consistent in the sense that for every 
state the SHF problem is solved independently by min- 
imization of the total energy. As seen from the figure, the 
TOMF impact is negligible. At the same time we have 
strong rearrangement effect for lp2h states, which is only 
accounted for the self-consistent mean-field calculations. 
This makes the simple estimations less reliable, although 
they still yield a nice first guess. 

Comparison with the experimental spectra shows that 
the calculated levels below 2 MeV follow the right se- 
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Fig. 15. The same as in Fig. [I4]but for the neutron spectrum 
in 133 Sn and the option min + s 2 alone. The dotted line marks 
the lowest collective state 2 + in reference even-even core 132 Sn. 
The experimental data are taken from [52] - 



quence and rather well reproduce the experimental ener- 
gies, especially for SLy6 and SV-bas. Discrepancies start 
for the 7/2~ state which has too high energy in all pre- 
dictions. Most probably this is caused by omitting the vi- 
brational admixture which generally lower and compress 
the low-energy spectrum. Figure Q3] shows the energy of 
the lowest collective state 3~ in the even-even core 208 Pb, 
which can serve as a rough energy indicator of the onset 
of core vibrational admixture. (A more complete estimate 
would, of course, require to know also the coupling matrix 
elements.) 

Results of similar calculations for neutron-excess nu- 
cleus 133 Sn are shown in Fig. [15j This is the case where 
the neutron particle (lp) states in 132 Sn are explored. Dif- 
ferent TOMF options are again found to make a negligible 
difference. Thus we display here only results from the final 
option min + s 2 . The simple lp estimate provides a nice 
first guess whose reliability, however, fades with increas- 
ing excitation energy, particularly if one approaches the 
vicinity of the 2 + core excitation. The comparison with 
the experimental energies in 133 Sn is less convincing than 
in 207 Pb. All excitation energies are overestimated. This 
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Fig. 16. The same as in Fig. [15] for the proton spectrum in 
131 In. The reference even-even core is 132 Sn. The experimental 
data are taken from 1531. 



is not a problem of spectral density because all three pa- 
rameterizations yield about the same trends. The most 
probable reason of the mismatch could be neglecting the 
core vibrational admixture whose impact in 133 Sn might 
be stronger than in 207 Pb and so affect even the lowest ex- 
citations. The comparison of the simple estimations with 
the min + s 2 case exhibits a strong rearrangement effect 
pertinent to self-consistent calculations. 

Fig. [16] shows the spectra for 131 In, the proton-hole 
neighbor of 132 Sn. The general trends are very similar to 
the previous figure. Again the simple lh spectrum pro- 
vides a good first guess, there is a strong rearrangement 
effect, and the calculated energies overestimate the exper- 
imental data. The later problem appears for all neighbors 
of 132 Sn while the neighbors of 208 Pb perform fairly well. 
This systematic deviation is still awaiting an explanation. 



4 Conclusions 

In this paper, we have discussed the description of odd- 
even and even-odd nuclei in the framework of the Skyrme- 
Hartree-Fock (SHF) method with BCS pairing. Particular 
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attention was paid to the influence of the time-odd mean 
fields (TOMF) in the Skyrme functional which are brought 
into action by the breaking of time-reversal symmetry 
through the odd nucleon. Three options were compared: 
even with time-even terms only and omitted TOMF, min 
with adding the minimum of TOMF to restore the Galilean 
invariance, and min + s 2 with additional spin couplings. 
Gradient spin couplings were omitted as leading to in- 
stabilities. The odd nuclei were computed within BCS 
with blocking the state occupied by the odd nucleon. A 
bunch of blocked states was computed self-consistently 
thus yielding the ground state and first few excitations in 
the odd nucleus. The nuclear properties which might be 
affected by TOMF, binding energies (BE), odd-even stag- 
gering (OES), separation energies (SE) and low-energy 
single-particle spectra were systematically investigated with 
three Skyrme parameterizations SkI3, SLy6, and SV-bas 
for the nuclear chart with 16 < Z < 92. 

The TOMF effects are found generally small and for 
the spectra even negligible. They affect the results much 
weaker than the choice of the Skyrme parameterization 
or rearrangement in self-consistent calculations. For BE, 
the TOMF effect is maximal for light nuclei and rapidly 
decreases with increasing nuclear mass number, in accor- 
dance to the previous results P^ll^6] . The change in BE is 
100 - 300 keV and TOMF may lead to both less and more 
binding, depending on the Skyrme parameterizations. The 
BE description is best with the minimal Galilean invariant 
option min while OES and SE prefer the full functional 
option min + s 2 . The OES can be modified by ~ 200 
keV and SE by ~ 200 keV. These results are fluctuations 
and sometimes conflicting since the currents and spin den- 
sities give opposite contributions and partly compensate 
each other, thus producing a fragile balance. Obviously, so 
weak and unstable TOMF effects cannot be used for the 
upgrade of the Skyrme functional in the spin channel. 

Some additional important points should be mentioned. 
First of all, the primary importance of the proper choice of 
the Skyrme parameterization has to be emphasized. This 
affects the results much more than TOMF. For all the nu- 
clear characteristics considered in the present exploration, 
the best description was obtained with the recent param- 
eterization SV-bas. Note that SV-bas has a large effective 
mass m*/m — 0.9 and its fit involves isotopic as well as 
isotonic chains. 

Second, differences of binding energies, like OES and 
SE, are especially sensitive to the shell and pairing de- 
tails. The OES is dominated by pairing in the mid shell 
regions and by shell effects near shell closures. Thus a 
strong dependence on the SHF parameterization is seen 
at shell closures and very little in mid shell regions where 
the agreement with experimental OES is generally satis- 
fying. It ought to be reminded that the pairing strength is 
calibrated to the spectral gap which deviates from OES. 
Hence the relation between these pairing characteristics 
was analyzed. 

Third, low lying excitation spectra for odd nuclei next 
to doubly magic ones have been analyzed. There are cases 
which work nicely and others which do not fit so well. The 



mismatch is probably caused by the coupling with vibra- 
tional and rotational modes of the even-even core. They 
are especially strong in deformed and transitional nuclei. 
The correlations could considerably improve the descrip- 
tion of BE, OES, and SE, first of all in mid shell nuclei. 
Besides, the correlations could be of a crucial importance 
for description of the low-energy single-particle spectra of 
odd nuclei. Being too dilute in the mean-field picture, the 
spectra can become more compressed and closer to the 
experimental data. 
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A Relation between parameters of the 
functional 



The SHF method can be alternatively formulated by start- 
ing since the density-dependent zero-range Skyrme force 

Vskyrme = Vl2 + V dd + V LS (21) 

V12 =io(l + x P CT )<5(r 1 -r 2 ) (22) 
+ |(l+x 1 P CT ) ( 5(r 1 -r 2 )fc 2 

+ t ±(\+ Xl P a )k 2 S{r 1 -r 2 ) 
+ t 2 (l + x 2 P a )kd(r 1 -r 2 )fe 
V& J = |(l + a:3P t r)p a ( EL |^)«(ri-r 2 ) (23) 



V LS = iti(<ri + <x 2 ) ■ k x <5(n - r 2 )k 



(24) 



with the spin-exchange operator P a = ^(1 + <Ti<t 2 ), the 
difference momentum k = — | ( Vi — V 2 ) acting to the 



right, and its counterpart k = | ^Vi — V 2 J acting to 

the left. Then the SHF functional ((HUH]) is obtained as the 
expectation value of this force with a Slater state. This 
yields one-to-one correspondence of the force parameters 
ti, Xi with the parameters b^b'^bi, b\ from ([6][9]): 



b 


= < (1 + hxo) 
= ±(l + \x 1 ) 

+ 1(1 + 1*2) 


b'o 


= t (h + 

- ^fi 4 
4 \2 \ 


Xo) 


h 


b[ 


Xl) 






*2 ( 1 

4 ^1 


f + x 2 ) 




= MKi + K> 

- 1(1 + 1*2) 


b' 2 


_ 3ti (1 
8 }2 


f SBl) 






+ tjL (- 
1 8 1 


r + X 2 ) 


h 


= 1*3(1 + 1*3) 
= 


b'. 


= |*3( 2 - 
= 3*4 


+ a; 3 ) 


h 


b', 




b 


= jhx 


b'o 


= 4*0 




h 


= §(*1&1 + t 2 x 2 ) 


b'i 


= I(*l " 


*2) 




= jq(MiXi - t 2 x 2 ) 
= 53*32:3 


b'2 


= 16(3*1 
= 2l*3 


+ < 2 ) 


h 


K 





(25) 
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B Single-particle wave function 

The single-particle wave function in the cylindrical coor- 
dinates is the spinor 

#<r)= (4-»wj = l4- , UU-"J' (6) 

with 

m[ s) = K k - ±s , m { - s) =K k + \s (27) 

and Kk is the projection of the total momentum of k-state 
on the symmetry z-axis. 
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